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ABSTRACT 



Aims. We model the X-ray reprocessing from a strong co-rotating flare above an accretion disk in active galactic nuclei. By performing 
detailed radiative transfer computations we explore the horizontal structure and evolution of the underlying hot spot. The main goal 
is to study how the resulting spectral features manifest themselves in short exposure time spectra. 

Methods. We analyze both the vertical and the horizontal spot structure and its dynamical reprocessed spectrum. To obtain the spectral 
evolution seen by a distant observer, we apply a general relativity ray-tracing technique. We concentrate on the energy band around 
the iron K-line, where the relativistic effects are most pronounced. Persistent flares lasting for a significant fraction of the orbital time 
scale and short, transient flares are considered. 

Results. In our time-resolved analysis, the spectra recorded by a distant observer depend on the position of the flare/spot with respect 
to the central black hole. If the flare duration significantly exceeds the light travel time across the spot, then the spot horizontal 
stratification is unimportant. On the other hand, if the flare duration is comparable to the light travel time across the spot radius, the 
lightcurves exhibit a typical asymmetry in their time profiles. The sequence of dynamical spectra proceeds from more strongly to 
less strongly ionized re-emission. At all locations within the spot the spectral intensity increases towards edge-on emission angles, 
revealing the limb brightening effect. 

Conclusions. Future X-ray observatories with significantly larger effective collecting areas will enable to spectroscopically map out 
the azimuthal irradiation structure of the accretion disk and to localize persistent flares. If the hot spot is not located too close to the 
marginally stable orbit of the black hole, it will be possible to probe the reflecting medium via the sub-structure of the iron K-line. 
Indications for transient flares will only be obtained from analyzing the observed lightcurves on the gravitational time scale of the 
accreting supermassive black hole. 

Key words, radiative transfer - accretion, accretion disks - galaxies: active - galaxies: Seyfert - X-rays: galaxies 
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The broad shape of the iron line complex seen in about 25% 
$H of all active galactic nuclei (AGN) is commonly assumed to 
5^ , be due to relativistic effects acting in the immediate vicinity 
of th e central, supermassive black hole ( e.g.'Tanaka et al."l995t 
Iwasawa et al. 1999; Wilms et al] 12001 : lLubiriski & Zdziarski 
2001; Fabian_^taLI I2002t lYoung et al. 1 120051: iGuainazzi et al.1 
2006). The line is expected to be a reprocessing feature emitted 
by the surface of the irradiated inner accretion disk. An attrac- 
tive scenario for the origin of the X-ray irradiation envisions the 
reconnection of magnetic loops above the disk surface (see e.g. 
Galeev, Rosner & Vayana 1979; Abramowicz et al. 1991; Haardt 
et al. 1994; Collin et al. 2003, Merloni & Fabian 2001). The hne 
emitting region lies at distances of only a few gravitational radii 
^ from the disk center. 
Recently, some spectrally resolved variabil ity of the line 
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comp l ex was observed fo r several obje c ts (i.e 
20041; jTurner et al.' '2006'; iMiller et al.l l2006t 
20071 Miniutti et al. 2007). For two of them, NGC 3516 and 
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NGC 3783, the iron line feature appears to vary systematically in 
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flux and centroid energy. The variations happen on characteristic 
time-scales of 25 ks a nd 27 ks respectively (I wasawa et al. 20041 
iTombesi et al.ll200'7l) . llwasawa et al.l (l2004l) suggest a promising 
explanation for this type of line behavior invoking an X-ray flare, 
which is co-rotating with the accretion disk. The flare irradiates 
the disk surface creating a hot spot. In the case of NGC 35 16 the 
flare and the underlying hot spot are located at (7-16) R^, from 
the black hole. The systematic changes of the reprocessed spec- 
trum are then explained by the orbital motion of the disk and the 
resulting time-dependent relativistic modifications. 

The above mentioned model requires a prescription for the 
reprocessed emission from the disk. The modeling of the re- 
processed spectra is often caiiied out with considerable sim- 
plifications. Many computations assume either a uniform den- 
sity of the disk atmosphere, or semi-isotropic illumination and 
"observation". The angular dependence of the reflected com- 
ponent formed by a constant density partially ionized medium 
was studied starting with Matt, Fabian & Ross (1993a,b), 
Zycki et al. (1994), Zycki & Czerny (1994). Then, reprocess- 
ing from a medium in hydrostatic equilibrium was modeled 
by Raymond (1993), Nayakshin et al. (2000), Ballantyne et 
al. (2001), Ballantyne & Ross (2002), Rozariska et al. (2002). 
Nayakshin & Kazanas (2002) pioneered the simplified time- 
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dependent reflection stu dies for photoionized accretion disks. 
iGoosmann et alj ([2006a') studied rms spectral variability from a 
distribution of flares with vertically stratified spots underneath. 
None of those models considered a horizontal stratification of 
the illuminated spot as suggested by Goosmann et alj (l2006bl) . 

In this paper we add considerable sophistication to the flare 
model by studying in detail the Compton reflection/reprocessed 
component coming from a single orbiting spot. Since the spot 
properties have a significant gradient from the spot center to 
the border, we take both the vertical stratification of the repro- 
cessing material and the horizontal stratification into account. 
The locally re-emitted spectrum is computed at various emis- 
sion angles, and combined with a rel ativistic ray-tracing method 
jDovciak et al.ll2004l; Id ovciaki 120041) to model the spectra seen 
by a distant observer. 

We model the observed time evolution of the hot spot emis- 
sion for two cases: a persisting flare lasting for a significant frac- 
tion of one orbital period at the given disk radius and short-term 
flares with durations of a few hundreds of seconds. The first case 
applies to recently reported observations of exceptionally power- 
ful flares lasting for minutes to hours. Example s were given e.g. 
for MCG-6-3 0-15 (see Fig. 7 in | Ponti et al.l2004 ) and NGC 5548 
(see Fig. 1 in lKaastra et al.ll2004l) . For the short-term flares, the 
horizontal stratification of the spot becomes important. During 
the onset and fading of the flare, the spot evolves from the center, 
where the first primary photons arrive, to the border. We calcu- 
late this effect presenting the resulting lightcurves and averaged 
spectra for short-term flares. We show that future data from the 
planned XEUS satellite will be of enough quality to be compared 
with our models. 

The structure of the paper is as follows; section |2] describes 
our model. The obtained locally emitted spectra are presented 
in Sect. |3] and those seen by a distant observer in Sect. |4] The 
importance of the studied effects is discussed in Sect.|5] Broader 
implications and some final conclusions are given in Sect. |6] 

2. The model 

Our method to compute the local spectra of the hot spot is simi- 
lar to the one described in Collin et al. (2003). First, we compute 
the vertical density structure of the disk at a given radius R as- 
suming that the disk is in hydrostatic equilibrium without any 
external irradiation. The density structure of the surface layers 
are then used to conduct detailed radiative transfer computations. 
We take into account the intrinsic time evolution of a flare and 
we compute the spectral changes as seen by a distant observer. 
We apply a ray-tracing technique to include appropriate time- 
delays, general relativity effects, and the orbital motion of the 
hot spot. 

In this paper, we assume for all computations the metric of 
a Schwarzschild black hole with M = lO^M© and an accretion 
rate ih = 0.001 in units of the Eddington accretion rate. The sur- 
rounding disk with the viscosity parameter a = 0.1 is considered 
at the radii 7 R„ and 18 R„. 

2. 1 . The vertical disk structure 

Following the argument made by Collin et al. (2003) we assume 
that the flare is much shorter than the characteristic timescale for 
the restoration of the hydrostatic equilibrium. Hence, we assume 
that there is no external irradiation before the onset of the flare 
and the computation of the initial disk profiles is based only on 
the emission coming from the standard a-disk at the given radius 



R. This condition is safely fulfilled if we assume that the flare 
duration is less than the orbital timescale at R. 

We assume that the disk atmosphere is a plane-parallel, verti- 
cally stratified slab. The density, temperature, and radiation flux 
profiles ar e computed using an a dvanc e d version of the code de - 
scribed in iRozaiiska & Czernvl (Il996l) . iRozariska et al. I (Il999l) . 
and Rdzariska et al. (2002). The code computes the hydrostatic 
balance of the medium in the gravitational field of the black 
hole, thereby including the relativisti c corrections to the equa- 
tion s of the vertical disk structure ( Novikov & Thorn"3 119731 : 
Pag e & Thornell974l) . The radiation inside the medium is trans- 
ported by difliision using the Rossland approximation. 

2.2. Computing fhe local spectra 

Having obtained the vertical density profile of the disk, we com- 
pute the local reprocessed spectra expected after the onset of the 
flare. Since the principal spectral features are produced within a 
few Thomson optical depths tes from the disk surface, we only 
consider the upper layers up to tes = 7 for the radiative transfer 
calculations. We call this upper layer the disk atmosphere and as- 
sume that it is heated from underneath by the thermal radiation 
of the underlying disk bulk. The temperature and incident flux 
^^disk of the disk bulk are computed when obtaining the vertical 
disk structure. From above, the disk atmosphere is irradiated by 
the flare located at a given height h that is measured from the 
disk surface where the flare causes the incident flux Fi„c- 

We assume that half of the primary radiation reaches the disk 
atmosphere. The main fraction of the flux is reprocessed directly 
underneath the flare source. The radiation emitted further side- 
ways reaches the disk farther away from the flare and so it be- 
comes geometrically diluted. Therefore, we can define a limited 
hot spot on the disk surface. The spot is centered on the orthog- 
onal projection of the point-like flare source. The circumference 
of the spot and the flare source define an irradiation cone with 
half-opening angle Oq. We neglect any irradiation outside this 
cone, thus we only consider reprocessed radiation coming from 
inside the spot. 

The distance R of the flare from the central black hole is mea- 
sured for the spot center We limit the flare cone to 6*0 = 60° and 
the height of the source to h = Q.SRg. We thus assume that the 
spot size is small and neglect the change in the disk's internal 
properties across the spot. It remains to be seen whether such an 
assumption is justified. The MHD simulations of Miller & Stone 
(2000) and Blaes, Hirose & Krolik (2007) suggest magnetic loop 
heights that are by ~ 2 - 3 times higher than the disk thickness, 
while Turner (2004) found little evidence for magnetic loops ris- 
ing so high. 

The total incident flux, Ff°J;, of the whole spot is parame- 
terized by the ratio between the total incident flux and the in- 
ternal flux emitted by the disk, Fdisk- We set Ff°J;/Fdi,sk = 144. 
This ratio ensures that the incident radiation by the flare is much 
stronger than the disk emission (the same value was used earlier 
by Ballantyne et al. 2001 and Collin et al. 2003). The local inci- 
dent flux, Fine, decreases with the distance from the spot center 
and the incident angle, 0i increases as illustrated in Fig. [1] We 
divide the hot spot into five concentric rings defining equally 
spaced solid angles with respect to the flare source. Integrated 
over its surface each ring thus receives the same energy per unit 
time. We solve the radiative transfer for each ring individually 
using separate values of Fine and di. Across a given ring we as- 
sume that Fine and 6i are constant. 

The spectral shape of the primary radiation is modeled as 
a power law extending from 1 eV to 100 keV with the pho- 
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Fig. 1. Illustration of the hot spot at the disk surface underlying 
the flare. The illuminated region is divided into concentric rings 
of equally intercepted fractions of the flare luminosity. 



ther away. Therefore there is a relative drift of the material across 
the spot, with some material leaving the spot region and cooling 
off' while new material enters and is heated. We assume that at 
each part of the spot the orbital motion is Keplerian but, for sim- 
plicity, we neglect the relative drift of the material and preserve 
the spot's circular shape We thus suppose that the heating and 
cooling of the drifting material is instantaneous. This is a safe 
assumption following the estimates by Collin et al. (2003). 

In this paper, we consider the relativistic light bending only 
for the radiation being re-emitted by the hot spot. Modeling the 
relativistic effects for the incident radiation is not necessary be- 
cause the flare occurs at a small height. We checked this ap- 
proximation conducting ray-tracing computations for incident 
photons being emitted along the limiting cone of the flare at 
9q - 60°. In flat space, these photons trace out a circle on the 
surface plane of the disk. For a flare in Schwarzschild metric at 
radius 7 /?g and height 0.5 the relativistic trace does not differ 
significantly from the flat-space result. 



ton index F = 1.9. The radiative transfer simulations are then 
conducted by the cou pled use of the codes Tit an and Noar 
described by Dumont, Abrassart. & CoUinl (l2000h and updated 
in IDumont et al.l (12003 '). The code Titan calculates the verti- 
cal ionization and temperature structure of the X-ray illumi- 
nated gas and the reprocessed spectrum in the optical, UV, 
and X-ray ranges. The radiative transfer equation is solved in 
multi-stream approximation by the Accelerated Lambda itera- 
tion (ALI) method for the continuum and for the lines. The 
resulting ionization and temperature data are then read by the 
Monte-Carlo code Noar. It calculates heating and cooling rates 
due to Compton scattering in the medium and the final emerg- 
ing X-ray spectra for the reprocessed radiation in the range of 
(0.8-100) keV. The resulting spectra are computed at 20 differ- 
ent emission angles ijj measured with respect to the disk normal 
and varying in constant steps of cos xf/. 



2.3. Keplerian motion and relativistic effects 

A distant observer detects an evolving spectrum to which all 
parts of the orbiting hot spot contribute. This contribution de- 
pends on the local spectral emissivity but also on the relativistic 
effects and on the Doppler shift acting at a given position. 

The radiation propagation for a Keplerian accretion disk 
arou nd a black hole is computed by the code KY (Dovci ak et al.l 
12004; Dovcia k2004). In our application, KY integrates the pho- 
ton geodesies between a position inside the orbiting hot spot and 
an observer at infinity. The integration method takes into account 
the energy shift of the photon along the geodesies, arrival time- 
lags, and the lensing-effect. Spectra are computed for various 
viewing directions, /, of a distant observer based on a given time- 
dependent emissivity distribution of the flare spot. We define the 
initial position of the spot center in polar coordinates Rq and 0o- 
With h = 0.5 Rf, and Go = 60°, the spot radius, ro, follows to be 

ro = V3/z = 0.866 

The flare above the disk surface is assumed to co-rotate with 
the underlying disk. This is a natural assumption if the flares 
form around reconnecting magnetic field lines that are anchored 
in the disk interior. The spot center follows the same Keplerian 
orbit. Since the spot originates from external illumination, we 
can assume that its initial shape remains circular. Nevertheless, 
in principle the kinematics of the spot is affected by the dif- 
ferential rotation of the disk: spot regions located closer to the 
black hole have a larger angular velocity than spot regions far- 



2.4. Time evolution of flare and spot 

The time evolution of the flare is based on two aspects. First, 
there is an intrinsic time evolution of the flare source, which 
is parameterized by a specific profile. For simplicity we as- 
sume that the intrinsic flare emission has a rectangular shape 
of length T. The source is thus switched on and off instanta- 
neously. Second, there is a delayed response from the various 
parts of the spot, as the primary radiation emitted by the elevated 
source reaches each concentric ring of the hot spot consecutively. 
During onset and fading of the flare, the hot spot evolves from 
the center, where the first primary photons arrive, to the border. 
This evolution adds time-dependent effects to the spectrum that 
are different from those caused by relativity and by the orbital 
spot motion. Let the irradiation at the spot center (Rq, ^q) start 
at time to and end there at to + T. Then the local emissivity s at 
the photon energy E, given disk position (R, (p), emission angle 
i/r, and time t is determined by: 



s(E,R,(f>, t//,t) = 



s.pJE, R, <p,ifr), to<t-ri<to + T, 



0, 



otherwise. 



(1) 



with rj - ( + - h)/c, and the spatial distance be- 
tween the spot center and the point (R, <p) is measured by r = 
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R^ + R^^- 2RRocos((f) - (po)- The local spot emission as com- 
puted from the radiative transfer simulations is denoted by Sspot- 
If the duration of the flare is much longer than the light travel 
time across the spot, the uradiation can be considered as a sta- 
tionary process and the onset and fading phases are negligible. 
Such a time-steady hot spot is described by the local disk emis- 
sivity 



s{E,R,(p,i//) = 



Sspot{E,R,(p,i//), r < ro, 

0, otherwise. 



(2) 



Comparison with the observational data requires relativis- 
tic ray-tracing and time integration of the local spot emissivity 
while taking into account the orbital motion. These computa- 
tions are consistently carried out by KY. 

3. The local properties of the reprocessed spectra 

The overall response of the irradiated disk to the incident flux 
is in agreement with known results: a hot layer forms in the up- 
per disk atmosphere, being at a temperature comparable to the 
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Compton temperature (Raymond 1993, Ko & Kallman 1994, 
Nayakshin et al. 2000, Ballantyne et al. 2001, Rozanska et al. 
2002). At larger optical depth a rapid transition to colder, less 
ionized layers follows. Nevertheless, as we show in the follow- 
ing, the horizontal stratification of the hot spot we introduced in 
our model brings new aspects to the structure of the medium and 
to the obtained reprocessed spectra. 

3.1. Temperature profiles and reprocessed spectra 

Figure |2] (top graph) shows the vertical temperature structure of 
the disk after the onset of the flare. We plot the temperature as 
a function of the Thomson optical depth and for the medium 
underneath three of the five concentric spot rings (inner, inter- 
mediate, and outer region of the spot). The profiles are plotted 
versus the Thomson optical depth tes, which increases toward 
the equatorial plane of the disk. The hot surface layer of the 
disk atmosphere and the abrupt transition to the colder medium 
are visible. At the spot center the temperatures are higher, the 
hot skin is thicker, and the less ionized layers lie deeper in the 
medium. For a given spot ring, the heated layer at 7 R„ is nar- 
rower than at 1 8 /?g which is determined by the different density 
profiles as computed before irradiation by the flare. 

The local reprocessed spectra emitted by different spot rings 
are shown in Fig. |2] (bottom graph). Going from top to bottom, 
the three panels present results for the inner, intermediate, and 
outer spot ring respectively. The three intensity spectra of each 
panel correspond to emission angles of^- 20° (solid line), ij/ - 
60° (long-dashed line), and tj/ - 80° (short-dashed line). At the 
lowest 4i, corresponding to a Seyfert-1 view, the Compton hump 
over (10-60) keV appears similar for all rings. At higher i/r, the 
normalization of the emitted spectra increases and the spectral 
slope softens from the spot center toward the border For a given 
ring, the spectral normalization rises with ijj constituting a limb- 
brightening effect of the hot spot. Such behavior had been first 
seen in models of irra diated constant density media (Fig. 7 of 
IZvcki & Czer"nvlll994 . 

The limb-brightening can be partly understood from the ver- 
tical temperature structure of the medium: The upper layers of 
the disk surface are very hot, therefore highly ionized, and opti- 
cally thin. A local observer looking at the medium at a normal 
direction tfr ~ 0° could see down to the deeper, less ionized lay- 
ers, because along this line-of-sight the temperature drops sig- 
nificantly within Tes - 1- However, a local observer looking 
at the same medium from a higher inclination tp ~ 60° would 
rather see the highly ionized surface layers dominated by elec- 
tron scattering and exhibiting less absorption. Along this line- 
of-sight the extinction within tes = 1 is dominated by scattering 
and thus more intensity is preserved than for the observer at low 
inclinations. 

We also investigate the dependencies of the reprocessed con- 
tinuum on if/ by considering spectral intensity ratios over narrow 
energy bands at 4 keV and at 28 keV. The results are plotted in 
Fig. [3] The plot shows that the spectra become harder toward 
lower emission angles. This behavior varies quantitatively with 
the location inside the spot. At the spot center the irradiating flux 
is strongest and the softest spectra appear This is true for all 
possible emission angles and can be understood by the relatively 
thicker high-temperature layer at the disk surface. It prevents a 
large fraction of the incident photons from reaching the deeper, 
colder medium, where they can be absorbed. Their chance to es- 
cape from the atmosphere after a few electron scattering events 
is higher than at the spot border, where the medium absorbs more 
efficiently. 



3.2. The iron fluorescence complex 

The spectral appearance of the iron K-line complex also varies 
with tfr and with the position inside the hot spot. Generally, 
several components of the line complex are visible: one at 
6.4 keV representing weakly ionized/neutral iron, the compo- 
nent of helium-like iron at 6.7 keV, and sometimes the hydrogen- 
like iron at 6.9 keV. A K/3-line at 7.1 keV is visible when the 
weakly ionized component at 6.4 keV is strong. For a given spot 
ring, the lines emitted at high t// indicate stronger ionization than 
at low iff. The effect is analogous to the limb-brightening dis- 
cussed above and can be again understood from the temperature 
and ionization structure of the medium. 

The iron K-lines and also the soft X-ray emission lines of 
the local spectra we compute are particularly strong with respect 
to the continuum. This is partly explained by the fact that we 
only consider the pure reprocessed radiation and we leave aside 
any dilution by the primary component. Another reason is that 
in our model the relatively cold disk matter at roughly 10^ K 
is irradiated by hard X-rays corresponding to a temperature of 
(10^-10^) K. Since we assume that the flare duration is sig- 
nificantly shorter than the dynamical timescale of the medium, 
the density structure of the disk atmosphere does not change. 
For longer fl are durations, the me dium must thermally expand, 
as shown in ICzernv & Goosmann (.2004.) and Rozanska et al. 
(2002). During the e xpansion the densi ty decreases leading to 
weaker line emission ('Collin et al."'20()3|). This effect can be im- 
portant for flares completing more than a single orbit, as the one 
analyzed in NGC 3516 bv llwasawa et al.1 (|2004 . 

4. Spectra seen by a distant observer 

To obtain the reprocessed spectra seen by a distant observer, we 
investigate the effects of the intrinsic time evolution of a flare 
combined with general relativity and with the orbital motion in 
the vicinity of the black hole. We apply the ray-tracing code KY 
with the local single-ring spectra described in Sect. |3] to define 
the emissivity of the spot. Note that in our modeling we focus 
on the re-emission of the spot and we do not take into account 
the primary emission. We thus assume that the flare source is 
shielded from the observer's direct view and that it manifests 
itself mainly by its reflection. 

4.1. Persistent spots at different orbital phases 

In a first step, we consider the case when the flare duration is 
much longer than the light crossing time of the hot spot. We 
thus neglect the onset and the decay phase of the flare when the 
spot is only partly illuminated. In the following we refer to this 
type of flares as persistent flares. To define the local emission we 
use the time-steady spot approximation (see Eq. |2]i. Beside the 
location of the flare Rq, (po, h, and the size of the spot Oo, ro given 
in Sect. 12.31 we assume a flare duration equal to 1/8 of the local 
orbital period. 

4.1.1. Spectral appearance 

In Fig.|4]we plot the spot spectra of a persistent flare occurring 
at different orbital phases 0o- The time-integrated spectra over 
(1-50) keV are shown for three disk inclinations with respect 
to a distant observer: / = 30°, / = 60°, and / - 85° (top, middle 
and bottom panels respectively). The impact of general relativity 
effects and of the spot motion becomes more important with in- 
creasing /. It affects the normalization of the reprocessed spectra 
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as well as the shape of specific features. At both distances from 
the black hole considered here, the relativistic modifications are 
relevant, but they are more significant at 7 Rg than at 18 Rg, as 
expected. The spacetime curvature and the Doppler shifts are 
stronger when being closer to the black hole. At! Rg, the iron K- 
line complex cannot be properly decomposed into components 
for any 0o and /. An exception occurs at / = 60°, where shortly 
"after" the passage "behind" the black hole some structure of the 
line appears. Still, the deformations remain quite strong also in 
this case. 

The spectral variations of a time-steady spot orbiting around 
a black hole are not only due to relativistic Doppler shifting and 
gravitational lensing. To an important extent they are also caused 
by time delay effects that are due to variations in the "length" of 
different photon geodesies reaching t he observer at a given time 
in his/her reference frame ( Dovciak et al.l 120071) . Gravitational 
lensing is especially important at edge-on viewing directions, 
when the spot, the black hole, and the observer are (nearly) lined 
up. At / - 85°, the maximum of the spectrum therefore occurs 
when the spot passes "behind" the black hole, at (po = 180° 
(top middle box in each panel). At intermediate inclinations the 
Doppler shifting is more dominant. In contrast to the Newtonian 
case, the maximum (minimum) relativistic Doppler shift does 
not occur when the spot is approaching (receding from) the ob- 
server at maximum velocity; it rather appears around ~ 225° 
{(p ~ 45°), which agrees with the results shown in Fig.|4l 

4.1 .2. Impact of the spot structure and of the 
angle-dependent local emission 

We investigate the impact of the horizontal spot structure and the 
angle dependent emission on the results. Therefore we compute 
time-integrated, relativistically convolved spectra for the follow- 
ing three cases of the local spot emission: 

• including the detailed horizontal spot-structure and the in- 
dependence of the local spot emission, 

• averaging over the horizontal spot structure but keeping the 
dependence on i/', 

• averaging over the horizontal spot structure and over t//. 

We use the same model setup as for Fig. |4] and we define the 
specific case of 0o = 225°, i.e. the spot is on the approaching 
side of the disk at the orbital phase where the strongest rela- 
tivistic modifications are expected. The spot then proceeds up to 
cj) = 270°. The inclination of the observer is set to / - 60°. This 
may correspond to a maximum type-1 viewing angle for very 
luminous AGN (quasars). 

We obtain the following results: the differences between a 
horizontally structured spot and a uniform one are negligible. 
The spot structure leaves a trace in the spectral evolution only 
during the rising and fading phases of the flare. These phases 
are very short compared to the total flare life-time. The time- 
integrated spectrum is therefore dominated by the reprocessed 
radiation coming from the fully illuminated spot, which evi- 
dently can be well approximated by uniform emission. Accurate 
computations of the angular emissivity dependence are slightly 
more important since they significantly affect the normalization 
of the reflected component. Nevertheless, the spectral shape is 
not strongly influenced: when the spectra are renormalized, the 
differences between the cases a) and c) are less than 7% at all 
energies. 

We conclude that for a spot seen at / - 60° and orbiting at 
R - 7 Rg it is safe to simplify the computations for persistent 



flares by averaging over the spot expansion. The same holds for 
lower inclinations, for other orbital phases than 0o = 225°, and 
for higher values of R because in all these cases the impact of the 
relativistic and Doppler effects is smaller than for our test case. 

4.1 .3. Detecting the orbital phase of a revolving spot 

In recent years, observations with XMM-Newton have started 
to reveal the radial irradiation structure of the innermost re- 
gions in AGN accreti on disks. For the Seyfert galaxy NGC 35 16 
llwasawa et al.l (|2004|) reported time-dependent spectral data that 
are in agreement with the modeling of an orbiting hot spot be- 
tween 7 Rg and \6 Rg. The data suggest a periodic signal in 
the iron line flux, which can be interpreted as a long-term flare 
completing several orbits. Nevertheless, the sensitivity of XMM- 
Newton is not sufficient to clearly resolve sub-orbital features in 
the spectral evolution of NGC 3516. Such investigations, which 
allow to map out the azimuthal irradiation structure of the inner 
accretion disk, will become possible with the much higher effec- 
tive collecting areas of future satellite observatories like XEUS 
and Constellation-X. 

To illustrate the spectral appearance of sub-orbital X-ray sig- 
nals we simulate the observation of an orbiting hot spot apply- 
ing our persistent flare model. We adopt the response matri- 
ces of the XEUS TES detector as it is available in April 2007 
and of the EPIC PN camera aboard XMM-Newton. The simula- 
tions for both satellites are conducted in xspec v. 11.3.2 using a 
modified version of the kylcr model dDovciak et al.ll2004 ). The 
model allows to implement our local spot spectra described in 
Sect. [3] after they are averaged over the horizontal spot expan- 
sion. According to Sect. 14.1.2! this simplification does not in- 
duce spurious effects. With kylcr we compute the relativisti- 
cally blurred, time-integrated spectrum of the spot that we use 
to produce the simulated XEUS data. We define a persisting spot 
lasting for 1/8 of an orbit at /? = 7 7?g on the receding (^o - 45°) 
and on the approaching {(po = 225°) side of the disk. We assume 
a Schwarzschild metric and a disk inclination of / = 30°. The 
simulated photon flux of 5 x 10"^ photons/(s cm^ keV) in the 
variable red wing of the iron line co mplex is adjusted to th e ob- 
servation of NGC 3516 analyzed in llwasawa et alj ( |2004|) . We 
assume an observation time T = 7300 s, which corresponds to 
the flare life time assuming a black hole of 10** solar masses. The 
fake data are rebinned using the tool grppha to ensure that each 
bin contains a minimum of 50 (200) counts for the XMM-Newton 
(XEUS) "simulation". 

In Fig. |5]we plot the produced spectra. For XEUS the flares 
on the two sides of the accretion disk can be distinguished from 
each other as the normalization and the position of the iron K- 
line differ significantly. Also in the XMM-Newton spectrum these 
differences are visible but the much higher noise level does not 
allow the same clear conclusion. Note that in these simulations 
we assume that the primary radiation does not dilute the reflec- 
tion. This holds true only if the object is strongly reflection- 
dominated or if the primary radiation has been precisely re- 
moved in the data analysis. In the latter case additional errors 
are expected, which makes it even less likely that XMM-Newton 
resolves such sub-orbital features. 

4.1 .4. Probing the local spot spectrum 

The computations of the local spectra we use include the ver- 
tical disk structure and the radiative transfer in such a strati- 
fied medium. The resulting reprocessed spectra reveal several 
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Fig. 5. Simulated observations by XMM-Newton (top) and XEUS 
(bottom) of a relativistically broadened iron K-line emitted by 
localized X-ray sources that orbit at 7 /?g around a 10*^ solar mass 
black hole. The two hot spots are at opposite positions, on the 
receding (0o - 45°, red) and the approaching ((^o = 225°, black) 
side of the accretion disk. 
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Fig. 6. Best fit for the simulated XMM-Newton (top) and XEUS 
(bottom) data of an approaching hot spot orbiting at 1 8 around 
a 10** solar mass black hole. The best fit models are indicated by 
solid lines: "power law + 2 kygline components" for the XMM- 
Newton data and "power law + 3 kygline components" for the 
XEUS simulations. 



iron line components, which are then blurred by the relativis- 
tic effects. The subtleties of the resulting spectrum cannot be 
well resolved in present-day X-ray data, and the question arises 
whether they will be visible with more advanced instruments in 
the future. We investigate to which extend the spectral structure 
in the iron line band can be observationally resolved by compar- 
ing simulated data from XMM-Newton and XEUS for a persistent 
flare. 

We consider a hot spot lasting for 1/8 of the orbital period at 
R - 18/?g corresponding to an integration time of 30700 s for 
a 10^ solar mass black hole. At this disk radius from the black 
hole, the relativistic effects are still relevant but not as violent as 
close to the marginally stable orbit. We therefore expect to still 
detect some structure of the iron line in the fake observations. 
The flare is supposed to start on the approaching side of the ac- 
cretion disk, at (/>() = 225°, and the system is seen at an inclina- 
tion of ; - 30°. The simulated observations are then obtained as 
described in Sect. 14.1.3] They are shown in Fig.|6] 

For the purpose of data analysis, it is a common contem- 
porary approach to relativistically blur one or several narrow 
Gaussians to approximate the local profile of the iron line, 
whereas the continuum is modeled by an underlying power law 
(which is invariant under the relativistic effects). We investigate 
how well the simulated data of a persistent spot can be repro- 
duced by such phenomenological models. We fit the simulated 
data consecutively adding more relativistically broadened line 



components to the underlying continuum power law. The lines 
are implemented by the kygline that is also available in xspec. 
The disk emission pattern is defined exactly as for the data sim- 
ulations, but the parameters of the power law as well as the line 
centroid, width, and normalization are allowed to vary in the 
fitting process. The fitting is restricted to the spectral range of 
(4.0-7.5) keV to avoid the soft X-ray features and the strong 
iron absorption edge. 

As expected the resulting phenomenological fits are much 
better determined for the more precise XEUS data than for the 
XMM-Newton simulations (Fig. |6j. For XEUS it is possible to 
find a phenomenological fit for three distinguished components 
of the Ka-line complex (Table [T]i: neutral/weakly ionized iron 
(~ 6.4 keV), helium-like iron (~ 6.7 keV), and hydrogen-like 
iron (~ 6.9 keV). The attempt to also identify an individual 
KjS-line around 7.1 keV fails, instead the hydrogen-like com- 
ponent seems particularly broad indicating a possible blend. All 
line centroids are slightly redshifted with respect to the labora- 
tory centroid energies, also the hydrogen-like Ka/neutral KJ3 line 
when interpreted as a blend. Visually the residuals imply that the 
fit misses some curvature in the spectrum red- and bluewards of 
the broadened iron line. The intrinsic spectral shape across the 
iron-line band is more complex than a power law continuum plus 
relativistically blurred Gaussian lines can account for. In the lo- 
cal reprocessing spectra, the various components of the iron K- 
line complex seem to stand on a "pedestal" between 6 keV and 
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Model component 


Parameter 


XEUS 


XMM-Newton 


XMM-Newton 


Power law 


1 




-u. i iz 


-u.uy I 




N 


1.42x10""^ 


1.91x10-°^ 


1.96x10"°^ 


KYGLINE #1 


E^^ [keV] 


O. Jol 




D.4jfS 




iT[eV] 


12.7 13 


20.243 


1 80.95 




A' 


6.06x10-* 


6.44x10-'*' 


1.09x10-°^ 


KYGLINE #2 


[keV] 


6.647 


6.659 


- 




iT[eV] 


Al.yj i 


nil*; 

U. 1 ID 






A' 


4.19x10-* 


4.35x10-°^ 




KYGLINE #3 


Ec [keV] 


6.983 








o-[eV] 


59.741 








N 


8.85x10-°' 








Red. X- 


352/250 


27/33 


39/36 



Table 1. Best-fit parameters for simulated XEUS and XMM- 
Newton observations of a hot spot at the phase of approaching 
motion. Power law models combined with up to three relativis- 
tically blurred iron lines are compared. 



7 keV (see Fig.|2|. This broad feature is mainly due to the drop 
in continuum normalization at the iron K absorption edges above 
7.1 keV. The edges should therefore be considered in the fitting. 
Redward of the 6.4 keV line, some impact of a Comptonized line 
shoulder might also be relevant. 

In the case of XMM-Newton, the number of required sub- 
components is much less well-defined (Fig. |6]l. Models with 
one or two Gaussian components both lead to acceptable fits 
(Table [1]), but including both, a neutral/weakly ionized and a 
helium-like Ka-component deliver the slightly better fit. For the 
model that only consists of one line the line width is very large 
as expected. In the two-line model only the neutral/weakly ion- 
ized iron line is broad and the helium-like component appears 
very naiTOW. Any further line around (6.9-7.1) keV remains un- 
detected. 

An important result of the modeling is, however, that repeat- 
ing the same data simulations for XMM-Newton shows a signif- 
icant uncertainty in the fitting results and their interpretations. 
The results for XEUS on the other hand are stable. This docu- 
ments that the sensitivity of XMM-Newton is not high enough to 
draw reliable conclusions about the sub-structure of the broad 
iron line shape coming from a persistent flare at 187?g as it is 
considered here. 

4.2. Transient spots at different orbital pliases 

In a second step we investigate the reprocessed spectra of flares 
with a life-time that is comparable to the light-crossing time of 
the spot. We name such spots as transient and we investigate 
their emission at defined orbital phases. The irradiating primary 
emission is switched on and off instantaneously. The emission 
lasts for one gravitational time scale, T,, = which compares 
to the spot's light-crossing time by Tg ~ 1.15 ^. During T„ the 
spot proceeds by the azimuthal angle A(p(R) - (Rg/Ry -^. At the 
two disk radii considered here this amounts to A<p(7 R„) ~ 3° 
and A0(18/?g) ~ 0.75°. Hence, during the flare irradiation the 
spot remains nearly at the same orbital phase. For the initial az- 
imuthal position 00 of the flare we consider the values <po - 0°, 
(f>o = 90°, 00 = 180°, and 0o = 270°. We use the same radiative 
transfer solutions for the local reprocessed spectra as in Sect. [3] 
We then run KY to compute sets of relativistically blurred, time- 
dependent spectra at the viewing angles / - 30°, / - 60°, and 
i = 85°. 



4.2.1. Spectral appearance 

In Fig.|2]we present the obtained integrated spectra for the ob- 
served flare periods and the coiTesponding lightcurves at the dif- 
ferent azimuthal positions. The observed flare durations are of 
the order of a few hundreds of seconds. The plots show that 
line broadening is important also for transient flares. But in con- 
trast to the persistent spots the spot velocity projected along the 
observer's line-of-sight does not change significantly during the 
flare period; this holds true even at high disk inclinations as the 
spot only moves by very small values of A(p. We therefore con- 
clude that the line broadening is caused to a large extent by gen- 
eral relativity. Due to the expansion of the spot, the photons ar- 
riving at the distant observer at a given time originate at different 
positions on the accretion disk. Hence, their gravitational red- 
shift, which only depends on the distance of the emission point 
from the black hole, is distributed over a certain range and the 
line appears to be broadened. At 18 /?g (right panels of Fig. |7]) 
the gravitational line broadening is weaker and the iron K-line 
complex is resolved into its sub-components at all orbital phases. 
When comparing the iron K-lines for the same 0o but at different 
inclinations, the effects of the vertical ionization and tempera- 
ture profile of the disk become visible (see Sect. [3]l and indicate 
that toward higher inclinations the medium appears to be more 
strongly ionized. 

Some aspects of the spectral shape are consistent with the 
results obtained in Sect. 14.11 for a persistent hot spot complet- 
ing a significant fraction of the orbit. The Doppler shift of the 
Ka line complex, for instance, coiTesponds to the orbital phase 
of the flare. It increases with higher inclination due to a larger 
velocity component projected along the line-of-sight. For a flare 
occurring "behind" the black hole, i.e. for 0o = 180°, the lensing 
effect strongly modifies the spectrum and increases the overall 
flux. At this phase the iron K-line is entirely smeared out when 
seen at high inclinations. 



4.2.2. Lightcurves 

In our modeling of transient flares we consider the different 
light-travel times between the elevated source and various posi- 
tions inside the hot spot. The effect can be seen in the lightcurves 
shown in Fig. |7] The lightcurves represent the integrated flux 
between 2 keV and 50 keV. The observed flare period is not 
the same for different values of / and 0o- Physically, it depends 
on the difference between the "longest" and the "shortest" pho- 
ton geodesies connecting the spot with the distant observer at a 
given time in the observational frame. The geodesies are a func- 
tion of the spacetime geometry and the photon flux they transfer 
depends on the time-dependent emission across the spot. 

At low inclinations the lightcurves have a characteristic and 
regular shape showing the rising and fading of the flare emission. 
The shape is not at all box-like as assumed for the incident radia- 
tion illuminating the spot. The geometrical expansion of the spot 
induces different time-delays for the reprocessed radiation and 
so the onset and the fading of the spot emission are smoothed 
out. A significant time asymmetry is seen in most cases, with 
a shallower rise and a more rapid decay. Going toward higher 
inclinations, the shape of the light curve changes and gener- 
ally becomes more symmetric - with an exception for the phase 
00 = 180°, when the spot passes "behind" the black hole. The 
closest symmetry is reached at ~ 60°. For even larger / the flux 
decreases rapidly and the lightcurves become much flatter. 



8 



Goosmann et al.: Iron lines from transient versus persistent hot spots 



4.2.3. Spectral time evolution 

To investigate the spectral evolution during a short flare we 
present in Fig. |8] snapshots of a time-dependent spot spectrum. 
For the particular model shown we consider a spot passing "be- 
hind" the black hole (^o - 180°) and a disk inclination of 30°, 
which is typical for type-1 AGN. 

The development of the flare emission across the hot spot is 
visible in the spectral sequence. The irradiation sets on at the 
spot center where it produces the highest temperature and ion- 
ization states. Therefore, the iron K-line complex shows stronger 
components of helium- and hydrogen-like iron for the rising 
branch of the lightcurve. During the flare period these compo- 
nents become weaker and the Ko' and K/3 emission of less ion- 
ized/neutral iron increase. The development of the line compo- 
nents is more clearly visible for a spot atR- 18 /?g. At 7 Rg, i.e. 
close to the marginally stable orbit, the relativistic broadening 
of the line makes it difficult to distinguish between its individ- 
ual components. Nevertheless, a development of the line from 
the onset to the decay of the flare can be seen. The spot evolves 
from a small emission site at its center to a circle/ring with outer 
radius ro = 0.866 Rg. Hence, the gradient in gravitational and 
Doppler shifting across the illuminated spot increases during the 
spot's life-time, so that the width and the centroid of the iron 
K-line evolve. 



5. Importance of the studied effects 

The spot model we present in this paper contains accurately 
computed local spectra for the reprocessed flare emission. We 
include the vertical stratification of the accretion disk and the 
geometrical details of the flare reprocessing before applying the 
relativistic modifications and computing the observed spectra. 
In previous work, these effects have not all been considered. We 
discuss in this section their impact on the results and we estimate 
to which extend it is necessary to include them into the model- 
ing. 

5. 1 . Local emission from fhe spot 

We calculated the local spot spectrum taking into account the 
vertical density stratification of the disk. Our local spectra ex- 
hibit more than one strong component of the iron line and a 
complex reprocessed continuum. The continuum component in- 
cludes a series of iron K absorption edges above 7.1 keV and a 
Comptonized shoulder on the red side of the Ka-line at 6.4 keV. 
These features give rise to a broad continuum "pedestal" carry- 
ing the emission lines, which is not completely smeared out by 
the relativistic effects even close to the marginally stable circular 
orbit of a Schwarzschild black hole. 

From simulating future XEUS observations of persistent 
flares in AGN at 18 7?g, we conclude that the sub-structure of 
the iron line complex can be investigated probing the reflecting 
medium. Using a detailed reflection model for the data analysis 
instead of phenomenological models has the advantage that the 
particular spectral characteristics are consistently accounted for 
The normalization of the emission lines with respect to the con- 
tinuum stays on physical grounds and can be connected to global 
parameters of the system (e.g. ionization parameter, disk ra- 
dius). For the case of Suzaku observations of AGN the Compton 
hump around 30 ke V can be fitted coherently with the broad iron 
K-line complex as bot h features are due to reprocessing (e.g. 
lGeorge&Fabianlll991h . 



Careful computations of the radial structure of a spot are not 
necessary if the time integration of the data is longer than the 
light travel time across the spot. In this case the spot spectrum 
can be approximated as a uniform (average) spectrum which 
simplifies the calculations. Since the flares considered in the cur- 
rent paper are located at a small height above the disk, the spot 
radii are small and this assumption is likely to be satisfied in 
most cases. The radial spot structure is going to be more impor- 
tant for flares located at large heights above the disk, but in this 
case the flare is a non-local phenomenon, and both the disk ra- 
dial structure as well as the light bending between the flare and 
the disk surface might be important. Such a model becomes es- 
sentially identical to the lamp-post geometry, which is discussed 
in Sect.llll 

5.2. General relativity and Doppler effects 

We have included photon energy shifts due to Doppler and gen- 
eral relativity effects, as well as light bending and mutual time 
delays that influence the observed signal. These effects have a 
strong impact on the observed spectral features and should be 
taken into account when modeling reprocessed radiation at dis- 
tances smaller than 20 Rg from the black hole. We find, that 
the relativistic effects do not only affect the observed shape 
of the iron K-line emission but also the underlying continuum 
"pedestal" that appears in the iron line band (see Sect. l5.lt . 

The relativistic effects are important for the localization of 
the hard X-ray emission. Short exposure time spectra, which 
cover only a fraction of the Keplerian timescale of the co- 
rotating flare, clearly exhibit the dependency of the observed 
spectral shape on the spot position. The observer's viewing angle 
is also important, as it determines the line-of-sight motion of the 
disk material, i.e. the Doppler shifts. At very high inclinations 
the gravitational lensing dominates the relativistic effects at the 
moment when the hot spot passes behind the black hole. This 
can be a relevant case for Compton-thin type-2 AGN. 

5.3. Flare model versus lamp-post model 

The flare model is a viable scenario for the primary X-ray 
emission in objects with relativistically broadened iron lines. 
Multiple flares are expected if the dynamo action leads to the oc- 
casional emergence of magnetic loops filled with plasma above 
the disk surface. Another popular model assumes that a single 
variable source is located near the symmetry axis; it could be ten- 
tatively identified with a shock at the jet base (Henri & Pelletier 
1991). The existence of a single "lamp post" implies that the 
loop formation is suppressed and a coUimated outflow is consid- 
ered instead. The two schemes differ rather deeply in the under- 
lying physics (Martocchia et al. 1996, 2000; Miniutti & Fabian 
2004). 

Both geometries currently have observational support. A sin- 
gle source at variable heights is able to explain the apparent lack 
of variability of the line expected to accompany the variations of 
the primary (Miniutti & Fabian 2004; Larsson et al. 2007; but see 
Zycki & Niedzwiecki 2006 for a critical analysis of the model). 
The multiple flare scenario is supported by detections of pow- 
erful irradiation occ asionally illuminating the disk at large radii 
dTurner et al.ll200^ . and in the case of many co-existing flares 
the scenario can also explain the apparent sup pression (Czerny 
et al. 2 004) and the relative lack of variability dGoosmann et al.l 
l2006al) of the iron K-line. 
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Obviously the two geometries share common properties, if 
a lamp-post source is off-axis (Dabrowski & Lasenby 2001), 
but differences should be detectable in their time-dependence. 
The light traveling times between the flare source and the re- 
processing disk are significantly larger for a lamp-post ge- 
ometry than for the flare geometry assumed in this paper 
Such delays can be constrained from cross-correlation analy- 
sis betweenjhejightcurves at different energy bands as shown 
in iGoosmann et alj ( l2007l) . Furthermore, we have shown in 
Sect. 14.1.31 how the more accurate spectroscopic observations 
expected from future satellites such as XEUS and Constellation- 
X will enable to map out the azimuthal details of the irradiated 
accretion disk - the azimuthal pattern of the disk re-emission de- 
pends on the radial position of the flare source, and is symmetric 
only when the flare is located on the disk axis, i.e. in the strict 
lamp post case. 



6. Discussion and conclusions 

We present exemplary results from very advanced time- 
dependent modeling of the reprocessed emission by a hot spot 
created underneath a magnetic flare. We attempt to test the sen- 
sitivity of these models to various simplifying assumptions. This 
gives the basis to determine a suitable grid of local solutions, 
which is needed to model an arbitrary flare with reasonable ac- 
curacy and efficiency, having in mind future X-ray missions like 
XEUS and Constellation-X. Our spot model includes the time- 
dependent effects that are caused by the delayed reprocessing 
from various parts of a spot. We have also included the horizon- 
tal spot structure and the angular dependence of the local emis- 
sion. The local emission is modeled taking into account the ver- 
tical disk structure and the observed time-dependent spot spec- 
trum contains all relativistic effects. We summarize the most rel- 
evant results we obtain: 

• The vertical temperature profile of the disk underneath 
the hot spot and the resulting local reprocessed spectra 
vary significantly across the spot surface. The ionized 
skin of the disk is thicker and the spectrum softer at the 
spot center than at the spot border. The different compo- 
nents of the iron K-line also follow this temperature gradient. 

• At all positions inside the spot, the local spectra show a 
limb-brightening effect and spectral softening with increas- 
ing emission angles ^. The iron line complex develops 
toward higher ionization states as ij/ rises. 

• The relativistically modified spectra seen by a distant 
observer depend significantly on the flare position. As 
expected, the spectral distortion of the iron K-line becomes 
more important at smaller disk radii and at higher disk 
inclinations. 

• For persistent flares the radial dependence of the local 
spectra across the hot spot can be safely averaged over, but 
the dependence of the local emission angle still has a small 
impact on the resulting spectrum seen by a distant observer 

• Current observational techniques are only at the threshold to 
determine the azimuthal phase of persistent hot spots. Their 
detection would require a perfect subtraction of an eventual 
primary component. 



• Transient spots show a slight spectral development reflect- 
ing the evolution of the local emission from the spot center 
toward the border. They reveal a large variety of shapes in 
X-ray lightcurves depending on the spot position and the ob- 
server's inclination. 

Especially the horizontal time-dependent re-emission of the 
transient spot has not been assessed in previous papers. It turns 
out that it has a significant impact on the observed lightcurve 
of the reprocessed radiation. The rising and the decaying phase 
have characteristic shapes and they do not look symmetric at 
disk inclinations smaller than 60°. This result will be testable by 
future observation techniques with a high-resolution and a high 
time-throughput. On the other hand, we obtain only minor dif- 
ferences in the spectral time-evolution between the rise and the 
decay of the spot emission. Even with XEUS or Constellation- 
X they should remain undetected as the maximum observation 
time is set by the flare duration. In our case we assume it to be 
1 Tg, equating only a few hundreds of seconds for a black hole 
with 10** Mq. 

The relevant time-scale of a transient flare as we assume it 
in this work scales with the black hole mass. One could there- 
fore imagine that the spectroscopic analysis of transient flares 
should become possible for higher black hole masses observed 
in some quasars. But unfortunately, high mass quasars tend to 
be more distant and therefore their count rate in the X-ray range 
is low compared to nearby, lower-mass Seyfert galaxies. Longer 
observable delay times between the primary and the reprocessed 
radiation are rather expected when the flare source is located 
higher above the disk than H - 0.5 as assumed in our model. 
But such a modeling approach requires to also include the rel- 
ativistic effects for the intrinsic radiation, which could be ne- 
glected for the small value of H we considered. 

The situation is much more promising for persistent flares. 
Using XEUS or Constellation-X to observe persistent spots in 
AGN the obtained data will be accurate enough to constrain 
not only the azimuthal location from X-ray spectroscopy (see 
Sect. 14.1.31 ), but it can even give some hints on the properties 
of the local spectrum. We find that the horizontal spot struc- 
ture has practically no impact on the observed spectrum so that 
the spot can be approximated by a uniform average emissivity. 
Nevertheless, this simplification would no longer be allowed for 
spots being close to the marginally stable orbit of a rotating black 
hole. In this case the relativistic modifications are more extreme 
and vary more strongly across the spot surface. 

In our current model, we have restricted ourselves to the case 
of flares that occur near above the disk (h ~ 0.5 GM/c^) and rel- 
atively far from the black hole (r > 6 GM/c^), so that the light- 
bending of the primary X-rays can be neglected. The primary 
reason for this restriction was motivated by our goal to concen- 
trate on the properties of the spot spectra in short-exposure ob- 
servations. This restriction means that we do not consider the 
influence of the black hole rotation, which would further com- 
plicate the analysis. Of course, a generalization to the rotating 
(Kerr) case should be explored in the future, because there are 
good chances that short-exposure spectra will remove the de- 
generacy which affects time-integrated spectral profiles. (There 
are no principal difficulties to include black hole rotation and our 
suite of numerical codes - Titan/Noar/KY - allows to treat the 
black hole angular momentum as a parameter.) 

For small distances to the black hole, light-bending effects 
must be important ( Miniutti & Fabian 2004) and the primary 
photons emerging from the flare source are directed down to 
the disk without reaching a very far distance from the flare. The 
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light-bending is an effect of the curved space-time in the direct 
vicinity of the black hole, and becomes weaker for flares at larger 
disk radii. We also assume that the intrinsic radiation is colli- 
mated inside an irradiation cone limited by = 60°. A possible 
interpretation of the irradiation cone can be given by assuming 
that the primary emission of the flare is not isotropic. A simple 
argument for such an assumption has been given by Nayakshin 
(2007) stating that magnetic flares occurring above the accre- 
tion disk of AGN are as likely to be anisotropic as solar flares. 
Nayakshin (2007) supports this argument with a simple consid- 
eration of the radiative transfer inside the coronal flare source. 
It is likely that at least for flares in the innermost region of the 
accretion disk both effects play a role. 

The results we present in this paper show that detailed radia- 
tive transfer calculations for the reprocessed flare emission and 
its modifications due to the relativistic effects in the vicinity of 
the black hole are computationally feasible with current technol- 
ogy. They deliver a vast variety of results for the time-dependent 
spectral energy distributions and for the flare lightcurves. To ap- 
ply them to a full extend to the observations of AGN, future 
generations of X-ray observatories are necessary. The models 
we presented here then allow to put important geometrical con- 
straints on the flaring activity in the accretion flow of supermas- 
sive black holes. 
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Fig. 2. The top graph shows the vertical temperature profiles for three concentric rings of the hot spot. They correspond to incident 
angles of = 18° (inner ring, top curves), 9 = 41° (intermediate ring, middle curves), and = 57° (outer ring, bottom curves). 
In the same order, the three panels of the bottom graph show the corresponding local reprocessed intensity spectra in £ x /g at 
three diff'erent emission angles: i// = 20° (blue, solid), i// - 60° (red, long dashes), and if/ = 80° (black, short dashes). The primary 
spectrum is indicated by the straight lines. The angles 6o and t// are measured with respect to the disk normal. Results for a hot spot 
slR^lRg (left) andR^lSRg (right) are presented. 
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Fig. 3. Spectral intensity ratios between 28 keV and 4 keV plotted versus the local emission angle ij/ (local disk inclination) for a hot 
spot at 7/?g (left) and at 18 7?g (right). The five curves denote the consecutive concentric rings of the hot spot. They are parameterized 
by the cosine of the incident angle 6 of the illuminating primary radiation. The angles 6 and i/' are measured with respect to the disk 
normal. 
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Fig. 4. Spectra emitted by a persistent hot spot at 8 different orbital phases. In each panel the bottom middle spectrum corresponds 
to the starting phase of the spot's closest approach to the observer (^o = 0°). The spot exists for 1/8 of the orbital time scale at its 
location and the disk rotates clock-wise. From top to bottom the panels represent different viewing angles, i (measured from the 
disk normal): i = 30°, i = 60°, and i = 85°. The spot is located at a distance of 7 Rg (left) and 18 /?g (right) from the black hole. 
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Fig. 7. Integrated flux spectra and lightcurves for short-term flares at four different initial azimuthal positions of the spot. From top 
to bottom, the panels represent 0o = 0°, (po - 90°, 00 = 180°, and 0o - 270° (continued on the next page). The origin of the 
azimuthal scale at = 0° is at the closest approach of the hot spot to the (inclined) observer; for (p - 180° the spot is farthest away 
and "behind" the black hole. The spectra were integrated over the entire flare period. The various curves of one panel denote the disk 
inclinations / - 30° (black, solid), / = 60°, (red, dashed), and / = 85° (blue, dot-dashed). The spectra contain a systematic offset 
by a factor of 5 between each other for clarity. Results for a spot at R - 1 Rg (left panels) and at R = 18 /?g (right panels) from the 
central black hole are shown. 
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Fig. 7. - continued 
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Fig. 8. Spectral evolution and lightcurves for a transient flare located "behind" the black hole (0o = 180°) with respect to the 
observer. The observer's inclination is / = 30°. The flare periods are divided into separate time intervals, over which the spectra are 
integrated. The corresponding time interval of each spectrum is marked on the flare lightcurve. Results for a spot at R - 1 Rg (left 
panels) and R - 18 (right panels) are shown. 



